Strategies for Optimize Off-Lattice Aggregate Simulations

We review some computer algorithms for the simulation of off-lattice clusters grown from a seed, with emphasis on the diffusion-limited aggregation, ballistic aggregation and Eden models. Only those methods which can be immediately extended to distinct off-lattice aggregation processes are discussed. The computer efficiencies of the distinct algorithms are compared.Comment: 6 pages, 7 figures and 3 tables; published at Brazilian Journal of Physics 38, march, 2008 (http://www.sbfisica.org.br/bjp/files/v38_81.pdf

On the relation between mass of pion, fundamental physical constants and cosmological parameters

In this article we reconsider the old mysterious relation, advocated by Dirac and Weinberg, between the mass of the pion, the fundamental physical constants, and the Hubble parameter. By introducing the cosmological density parameters, we show how the corresponding equation may be written in a form that is invariant with respect to the expansion of the Universe and without invoking a varying gravitational "constant", as was originaly proposed by Dirac. It is suggest that, through this relation, Nature gives a hint that virtual pions dominante the "content" of the quantum vacuum

Contact process on a Voronoi triangulation

We study the continuous absorbing-state phase transition in the contact process on the Voronoi-Delaunay lattice. The Voronoi construction is a natural way to introduce quenched coordination disorder in lattice models. We simulate the disordered system using the quasistationary simulation method and determine its critical exponents and moment ratios. Our results suggest that the critical behavior of the disordered system is unchanged with respect to that on a regular lattice, i.e., that of directed percolation

Aggregation in a mixture of Brownian and ballistic wandering particles

In this paper, we analyze the scaling properties of a model that has as limiting cases the diffusion-limited aggregation (DLA) and the ballistic aggregation (BA) models. This model allows us to control the radial and angular scaling of the patterns, as well as, their gap distributions. The particles added to the cluster can follow either ballistic trajectories, with probability $P_{ba}$, or random ones, with probability $P_{rw}=1-P_{ba}$. The patterns were characterized through several quantities, including those related to the radial and angular scaling. The fractal dimension as a function of $P_{ba}$ continuously increases from $d_f\approx 1.72$ (DLA dimensionality) for $P_{ba}=0$ to $d_f\approx 2$ (BA dimensionality) for $P_{ba}=1$. However, the lacunarity and the active zone width exhibt a distinct behavior: they are convex functions of $P_{ba}$ with a maximum at $P_{ba}\approx1/2$. Through the analysis of the angular correlation function, we found that the difference between the radial and angular exponents decreases continuously with increasing $P_{ba}$ and rapidly vanishes for $P_{ba}>1/2$, in agreement with recent results concerning the asymptotic scaling of DLA clusters.Comment: 7 pages, 6 figures. accepted for publication on PR

Is it really possible to grow isotropic on-lattice diffusion-limited aggregates?

In a recent paper (Bogoyavlenskiy V A 2002 \JPA \textbf{35} 2533), an algorithm aiming to generate isotropic clusters of the on-lattice diffusion-limited aggregation (DLA) model was proposed. The procedure consists of aggregation probabilities proportional to the squared number of occupied sites ($k^2$). In the present work, we analyzed this algorithm using the noise reduced version of the DLA model and large scale simulations. In the noiseless limit, instead of isotropic patterns, a $45^\circ$ ($30^\circ$) rotation in the anisotropy directions of the clusters grown on square (triangular) lattices was observed. A generalized algorithm, in which the aggregation probability is proportional to $k^\nu$, was proposed. The exponent $\nu$ has a nonuniversal critical value $\nu_c$, for which the patterns generated in the noiseless limit exhibit the original (axial) anisotropy for $\nu<\nu_c$ and the rotated one (diagonal) for $\nu>\nu_c$. The values $\nu_c = 1.395\pm0.005$ and $\nu_c = 0.82\pm 0.01$ were found for square and triangular lattices, respectively. Moreover, large scale simulations show that there are a nontrivial relation between noise reduction and anisotropy direction. The case $\nu=2$ (\bogo's rule) is an example where the patterns exhibit the axial anisotropy for small and the diagonal one for large noise reduction.Comment: 12 pages, 8 figure

Hawking Radiation in the Dilaton Gravity with a Non-Minimally Coupled Scalar Field

We discuss the two-dimensional dilaton gravity with a scalar field as the source matter where the coupling with the gravity is given, besides the minimal one, through an external field. This coupling generalizes the conformal anomaly in the same way as those found in recent literature, but with a diferent motivation. The modification to the Hawking radiation is calculated explicity and shows an additional term that introduces a dependence on the (effective) mass of the black-hole.Comment: 13 pages, latex file, no figures, to be published in IJM

Morphological transition between diffusion-limited and ballistic aggregation growth patterns

In this work, the transition between diffusion-limited and ballistic aggregation models was revisited using a model in which biased random walks simulate the particle trajectories. The bias is controlled by a parameter $\lambda$, which assumes the value $\lambda=0$ (1) for ballistic (diffusion-limited) aggregation model. Patterns growing from a single seed were considered. In order to simulate large clusters, a new efficient algorithm was developed. For $\lambda \ne 0$, the patterns are fractal on the small length scales, but homogeneous on the large ones. We evaluated the mean density of particles $\bar{\rho}$ in the region defined by a circle of radius $r$ centered at the initial seed. As a function of $r$, $\bar{\rho}$ reaches the asymptotic value $\rho_0(\lambda)$ following a power law $\bar{\rho}=\rho_0+Ar^{-\gamma}$ with a universal exponent $\gamma=0.46(2)$, independent of $\lambda$. The asymptotic value has the behavior $\rho_0\sim|1-\lambda|^\beta$, where $\beta= 0.26(1)$. The characteristic crossover length that determines the transition from DLA- to BA-like scaling regimes is given by $\xi\sim|1-\lambda|^{-\nu}$, where $\nu=0.61(1)$, while the cluster mass at the crossover follows a power law $M_\xi\sim|1 -\lambda|^{-\alpha}$, where $\alpha=0.97(2)$. We deduce the scaling relations \beta=\n u\gamma and $\beta=2\nu-\alpha$ between these exponents.Comment: 7 pages, 8 figure

Modelling of epitaxial film growth with a Ehrlich-Schwoebel barrier dependent on the step height

The formation of mounded surfaces in epitaxial growth is attributed to the presence of barriers against interlayer diffusion in the terrace edges, known as Ehrlich-Schwoebel (ES) barriers. We investigate a model for epitaxial growth using a ES barrier explicitly dependent on the step height. Our model has an intrinsic topological step barrier even in the absence of an explicit ES barrier. We show that mounded morphologies can be obtained even for a small barrier while a self-affine growth, consistent with the Villain-Lai-Das Sarma equation, is observed in absence of an explicit step barrier. The mounded surfaces are described by a super-roughness dynamical scaling characterized by locally smooth (faceted) surfaces and a global roughness exponent $\alpha>1$. The thin film limit is featured by surfaces with self-assembled three-dimensional structures having an aspect ratio (height/width) that may increase or decrease with temperature depending on the strength of step barrier.Comment: To appear in J. Phys. Cond. Matter; 3 movies as supplementary materia